A Molecular Platinum Cluster Junction: A Single-Molecule Switch 
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Abstract: We present a theoretical study of the electronic transport 
through single-molecule junctions incorporating a Pt6 metal cluster 
bound within an organic framework. We show that the insertion of 
this molecule between a pair of electrodes leads to a fully atomically- 
engineered nano-metallic device with high conductance at the Fermi 
level and two sequential high on/off switching states. The origin of 
this property can be traced back to the existence of a HOMO which 
consists of two degenerate and asymmetric orbitals, lying close in 
energy to the Fermi level of the metallic leads. Their degeneracy is 
broken when the molecule is contacted to the leads, giving rise to 
two resonances which become pinned close to the Fermi level and 
display destructive interference. 



Molecular metal clusters have recently started being ex- 
plored in breakjunction type single-molecule-conductance experi- 
ments, 1-3 whilst previous studies have also investigated clusters in 
the STM-STS configuration which mantains a tunnel gap between 
the tip and the molecule.^ Clusters have fewer metal atoms than 
metal nanoparticles but a more precisely defined composition and 
structure and as such they are better suited for comparative studies 
between theory and experiments. Bare metal clusters have been 
successfully deposited and studied on various surfaces. 6,7 Such 
clusters however require very clean conditions for study due to their 
high reactivity. In order to overcome this fact, and study molecules 
under the normal ambient conditions of the breakjunction experi- 
ment, the metal atoms must be encapsulated in a ligand sphere to 
avoid unwanted chemical reactions with the atmosphere. Ligand 
stabilized metal clusters (in particular those incorporated in an or- 
ganic framework) have therefore been suggested as components of 
data storage devices where they would act as nano-capacitors due 
to their redox properties. 8 Transport measurements on individual 
molecular metal clusters are, however, still rare. Metal atoms have 
been successfully incorporated into organic frameworks as metal 
complexes for single-molecule experiments, as well as chains of 
atoms. 9-12 Regarding metal clusters, the Mni2 structure has been 
investigated as the functional core in single molecule magnets 1 and 
theoretical studies about its transport properties have been carried 
out. 13,14 However, despite the vast range of molecular clusters 
known, only a few have been analyzed in this context. This is 
especially true regarding the theory of the electronic transport. 15 i 16 

In ref., 2 single-molecule experiments on 
[Pt 6 (M-P t Bu 2 ) 4 (CO) 4 (S(CH 2 ) 4 SH) 2 ) 3 ] (from now on (SC 4 S) 2 Pt 6 ) 
were carried out using the STM based I(s) technique. 17,18 The 
structure of this Pt cluster makes it ideal for these kind of measure- 
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ments as the ligands which act as binding groups (the alkanethiol 
chains) are oriented trans to each other, producing a linear wire. 
Without these chains, it would be difficult to determine a precise 
point of attachment between the molecule and the electrodes. The 
presence of these groups not only provides robust thiol anchors, 
but also increases the total length of the molecule in one defined 
direction, allowing measured break-off distances to assess accu- 
rately the orientation of the molecule in the junction (the need of 
incorporating the metal cluster in a wire has also been suggested in 
a theoretical study of the electronic transport through a Mni2 clus- 
ter 14 ). In the experimental study, it was suggested that the presence 
of the platinum cluster increases the conductance compared to an 
alkane chain of equivalent length due to there being Pt states in 
close proximity to the Fermi level. These were claimed to create 
an indentation in the potential barrier, giving a molecular analog 
of an inorganic double tunneling barrier, with the Ptg unit acting 
as a well and the alkane chains acting as barriers so to decouple 
the central unit from the leads. This was, however, a speculative 
argument which was never fully proven, and neither was the nature 
of these states identified. Finding out whether this assertion could 
be true is important because the molecular double tunneling barrier 
configuration could be the basis for creating molecular switches, 19 
analogous of conventional electronic components. This is ulti- 
mately one of the main aims of molecular electronics. 20 In this 
study presented here, we carried out first principle calculations 
to find out if we could corroborate the assumptions made in the 
experimental work. We will show that the HOMO of (SC 4 S) 2 Pt 6 
consists of two degenerate levels which get pinned to the Fermi 
level when the molecule is placed between two gold electrodes. 
However, these two levels do not originate simply from the Ptg 
unit, but contain states localized on the ligands also. Moreover, 
interesting interference features appear in the transmission curve at 
the Fermi level. 

We optimized all geometries with Turbomole 6.1, 21 using the 
BP86 functional 22 and the def-SVP basis setM We built the molec- 
ular junctions by placing the molecules (previously relaxed in the 
gas phase) between two 20 gold atom clusters and then relaxed the 
molecule and the 4(3) innermost gold atoms on each side in a top 
(hollow) geometry. The subsequent transport calculations were car- 
ried out with ANT, 24 which is built as an interface to Gaussian. 25 
This code computes the electronic transmission using nonequilib- 
rium Green's function techniques in the spirit of the Landau er for- 
malism, employing parametrized tight-binding Bethe lattices in the 
electrode description. 26 For these calculations, we used the PBE 
functional, 27 a lanl2DZ basis set for Pt atoms and the four inner- 
most Au atoms on each side, 28 while we used a CRENBS basis 
setS for all the other Au atoms. A 6-31++G basis set was chosen 
for C, S, P, O, and sto-3g for H. 

We first studied the molecule in the gas phase. The optimized 
geometry is shown in Q] It consists of a Ptg cluster sandwiched 
between two butanedithiol (four carbon atoms) chains. The cen- 
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Figure 2. (Color online) Optimized geometries for (SC4S)2Pt6 and Cig 
embedded between Au clusters. 



Figure 1. Central moiety of (SC4S)2Pt6 (a) and optimized geometry of the 
whole molecule with the two degenerate HOMO orbitals (b-c). 



tral cluster consists in turn of two orthogonal Pt3 triangles, con- 
firmed by crystallographic structural determination on the carbonyl 
substituted compound. 30 Each Pt3 triangle contains two bridging 
phosphine groups and the remaining sides are joined together to 
form a Pt tetrahedral cored. The tetrahedral core also contains four 
CO ligands. The molecule is slightly bent, with the angle between 
the alkyl chains less than 180°. In the second and third panel of 
[T]we show the HOMO. It consists of two doubly occupied levels 
which are degenerate and the corresponding orbitals of which are 
related to each other by a 4-fold improper rotation. It is mainly 
localized on the two innermost S atoms and on the two neighbor- 
ing Pt atoms. Notice that these frontier orbitals are different from 
the HOMO shown in ref 30 due to the presence of the butanedithiol 
ligands. We also studied a C18H36 alkyl chain terminated with S 
atoms (from now on C\g), as it was claimed 2 to have the same 
length as (SC4S)2Pt6 but a much lower conductance. Such a com- 
parison would help us decide whether the presence of the Pt cluster 
provides a significant change in the conductance with respect to the 
alkyl chain alone. We first checked that the length of this chain pro- 
vides a reasonable comparison to the cluster molecule. The length 
of C is is 24.9 A, which is halfway between the through space S to 
S distance (23.9 A) and the length measured by adding the 3 com- 
ponents (2 chains + Pt core, 25.8 A ) in (SC 4 S) 2 Pt 6 . In Cig, the 
HOMO is localized on the S atoms, while the HOMO-1 is delocal- 
ized throughout the whole chain— . 

After analyzing the molecules in the gas phase, we proceeded to 
study junctions containing (SC4S)2Pt6 and Cis bound to the gold 
clusters in a top and a hollow geometry (0. [3] shows the transmis- 
sion curves for all four cases. In the case of Cis, the energy align- 
ment seems to be considerably affected by the binding geometry. In 
the hollow geometry, the HOMO-1 (which is delocalized through- 
out the whole chain and gives the main contribution to the transmis- 
sion 20 ) is shifted down in energy compared to the top position (as 
has also been observed for other thiolated molecules 31 ), yielding a 



difference in the conductance of one order of magnitude between 
the two geometries. Notice also the presence of the HOMO (local- 
ized on the S atoms) very close to the Fermi level, clearly visible as 
a bump in the transmission curve for the top geometry. In the case 
of (SC4S)2Pt6, the alignment of the HOMO states does not show a 
particular dependence on the binding geometry as they are pinned 
at the Fermi level in both cases, but rather their splitting seems to 
be affected. The degeneracy of the HOMO states is broken in the 
junction as the symmetry has changed upon geometry relaxation. 
We stress that the two peaks close to the Fermi level do not arise 
simply from the Ptg unit as suggested in ref., 2 but rather from or- 
bitals which contain both S and Pt contribution, which however do 
not spread over all Pt atoms, as discussed above. 

Fermi level pinning has been widely studied but its nature is 
still under discussion (see ref— and references therein). It is also 
well known that molecular HOMO-LUMO gaps are usually un- 
derestimated in DFT and that this affects the conductance val- 
ues in DFT-based transmission calculations. In order to gain in- 
sight into the level alignment in our calculations and to understand 
whether it is reliable, we evaluated the ionization potential (IP) and 
electron affinity (EA) in the gas phase for the two molecules, as 
E(N) -E{N- 1) and E(N + 1) -E(N), respectively, where N is 
the total electron number. These quantities (based on total-energy 
differences) are expected to be more reliable than the HOMO and 
LUMO assigned by the Kohn-Sham states (at -4.27 and -3.07 eV 
for (SC 4 S) 2 Pt 6 and at -5.3 and -0.08 eV for Cig, respectively). The 
calculated IP and EA in the gas phase were found to be at -5.66 
and -1.70 eV for (SC 4 S) 2 Pt 6 , and -7.4 and 0.45 eV for Cig, re- 
spectively. Molecular HOMO-LUMO gaps are expected to narrow 
when molecules approach metal electrodes due to the screening ef- 
fect.— In the case of (SC4S)2Pt6, since the IP is very close to the Au 
Fermi level (which is at -5.0 eV), it is plausible that contacting the 
molecule to the electrode raises the IP until such a point that it be- 
comes pinned to the Fermi level. This is accompanied by a charge 
transfer from the molecule onto the metal, confirmed by 1.47 and 
4.3 e positive charge found on the molecule in the top and hollow 
position, respectively. Notice that the transferred charge originates 
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Figure 3. Transmission as a function of energy for all studied molecular 
junctions. 



mostly from the S atoms directly connected to the leads. Based 
on the proximity of the IP to the gold Fermi level, we believe that 
the energy alignment of the HOMO states appearing in our trans- 
mission curves for (SC4S)2Pt6 is robust and unaffected by common 
DFT failures, providing a realistic picture of the experimental sce- 
nario. On the other hand, in the case of Cis, the IP calculated in the 
gas phase is at a much lower energy than the peak corresponding to 
the HOMO (localized on the S atoms) in the transmission curve for 
the Au-Ci8-Au junction, especially for the top binding geometry. 
This suggests that it should align at lower energy; if so, the conduc- 
tance values of Ci8 would be significantly lower than those which 
we computed. 

A further significant result is the presence of interference res- 
onances in the transmission curves, which are clearly visible in [3] 
The interest in interference effects in molecular junctions is steadily 
growing 34-37 and this stems from the fact that they modify the ther- 
moelectric properties of the junctions 20 and that the ensuing dips in 
the transmission curves could give rise to large on/off ratios in fu- 
ture molecule-based electronic devices. 36 Interference effects have 



recently been detected experimentally in molecular junctions.— In 
our case, they originate from the existence of two possible pathways 
(the two HOMO states). For the top position, we find a similar fea- 
ture to that produced in ref. 39 Interferences yielding peaks at the 
Fermi level have already been predicted, 40,41 however their occur- 
rence from two degenerate levels pinned at the Fermi level, to the 
best of our knowledge, is rare. In order to shed light on the origin of 
these interference features, we built a simple model as depicted in[4] 
where we consider two levels both connected to the leads. They are 
related to the two HOMO levels of (SC4S)2Pt6, which are degener- 
ate and orthogonal in the gas phase. However, when the molecule 
is placed between two electrodes, the symmetry and consequently 
the degeneracy is broken due to the geometrical readjustment of 
the molecule. Thus, in our model, we consider two states which are 
no longer eigenstates of the system and are coupled by a non-zero 
hopping matrix element t. 

In this model, the Hamiltonian of the central (molecular) part 
takes the following form: 
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Figure 4. (a) Schematic representation (black) of the toy model describing 
the interference effects close to the Fermi energy and of the double tunnel 
barrier (blue), (b) The corresponding transmission as a function of energy 
for the case of the top position (left panel) with the following parameter 
values: F = 0.05 eV, e = -0.7 eV, £i = £ 2 = -0.12 eV, h = 0.009 eV, 
t 2 = 0.005 eV, t = 0.14 eV, tu = 0.003 eV, and for the hollow geometry 
(right panel) with the parameter values: F = 0.05 eV, £o = —0.7 eV, £\ = 
e 2 = -0.06 eV, h = 0.005 eV, t 2 = 0.002 eV, t = 0.03 eV, tu = 0.000145 
eV. 



coupled to each other via t. £qis the energy of the interface levels, 
coupled to £\ and £2 through t\ and *2- Notice that t\ must be bigger 
than t2, in order to account for the asymmetry of each orbital. The 
transmission is given by 

--*r L r R \G u \ 2 



T(E) 



(1) 



where 



where £\ and £2 are the energies of two degenerate molecular levels, 



G=[E-H C -L]- 1 (2) 

and E = El + Er, being El and Er the self energies for left and 
right lead, respectively. These two matrices have all terms equal to 
but Em = /Jl and ER44 = being II and the coupling to 
left and right lead, respectively. In our case, II = = F as we are 
considering symmetric junctions. 

As an example, in the inset of |4]we show the transmission curve 
given by this model for a chosen set of parameters (listed in the cap- 
tion) which can reproduce the shape of the DFT-calculated trans- 
mission curve for both hollow and top position. This proves that 
the nature of the features appearing in the transmission curves at 
the Fermi level is in the interference between the two levels of the 
central moiety and two levels localized at the interface between this 
central part and the leads. In the case of asymmetric junctions, such 
as top-hollow, the transmission curve is expected to resemble that 
corresponding to the top-top case (see the Supporting Information). 

As mentioned above, the presence of a resonance immediately 
followed by an antiresonance increases the on/off ratio in gating 
experiments. One can simply view this by imagining that the Fermi 
level (located at eV) is shifted to more negative values along the 
energy axis: by doing this, in the top geometry, for instance, the 
Fermi level will cross a first sharp peak (on) and, after a very small 
shift, a sharp dip (off). This sequence is then repeated by shift- 
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ing the Fermi level further so as to cross the second resonance- 
antiresonance couple. Indeed, we propose that the low bias con- 
ductance of this system should be measured experimentally as a 
function of the gate voltage, in order to test the breaking of the de- 
generacy. 

Finally, we compare our results with the experiments. In ref., 2 
a conductance of 3 x 10 -5 Go was measured for (SC4S)2Pt6, while 
the conductance of C\% was extrapolated from the experimental at- 
tenuation (beta) value measured for a series of shorter compounds 
as 6 x 10 _10 G . It was argued that the alkyl chains in (SC 4 S) 2 Pt 6 
simply act as spacers as the frontier orbitals lie far from the Fermi 
level, while the frontier orbitals of the central unit create a barrier 
indentation and, consequently, raise the conductance by approxi- 
mately 6 orders of magnitude compared to Cis- In our results, it 
is true that precisely at the very Fermi level the conductance of 
(SC4S)2Pt6 is definitely higher than Cis; however, at other ener- 
gies immediately below or above it, the conductance is compara- 
ble or even lower, especially as a result of the destructive interfer- 
ence. This could potentially explain some non-linearities observed 
in some of the I-V curves recorded experimentally (see the Sup- 
porting Information of ref—). It is also worth adding that, due to 
the uncertainty about the precise length of the Pt molecule, exactly 
comparing the conductances of (SC4S)2Pt6 and Cis is somehow 
ambiguous, since, in our equilibrium geometries, Cis is straight 
(with only slight defects) whereas the cluster molecule is bent, with 
the Ptg cluster slightly out of the Au-Au axis. In the experiments 
though, the Pt molecule is probably straighten, due to the pulling 
stress applied. Regardless these uncertainties, the overall picture 
supports the naive idea based on the experimental finding, confirm- 
ing the presence of resonances at the Fermi level arising from the 
presence of the central additional moiety in the molecular cluster. 

In conclusion, we have theoretically studied single molecule 
junctions incorporating a Ptg cluster, showing that it acts in an anal- 
ogous fashion to a double tunneling barrier due to two (originally 
degenerate) states which align at the Fermi level. These states do 
not, however, stem from the Pt unit alone, but specifically from two 
apical Pt atoms and their neighboring S atoms. This gives rise to 
quantum interference effects due to multiple electronic pathways 
through the molecule. Due to the pinning, these effects should be 
easily detectable experimentally. Despite the seemingly apparent 
complexity of the structure, all the physical properties arise from 
these four atoms alone, highlighting the delicate relationship be- 
tween molecular structure and electrical properties. We propose 
that further chemical synthesis combined with theoretical calcula- 
tions should be performed in order to explore more combinations of 
metal clusters with organic moieties. This could provide an alter- 
nate design strategy for molecular devices other than through purely 
organic molecules due to the wealth of different structures possible 
via inorganic chemistry. Recent progress in implementing a gate 
in single-molecule based devicesi 2 ^^ allows us to predict for the 
studied molecule, with its controlled binding to the electrodes and 
with its two levels aligned at the Fermi level, that it could be used in 
a three-terminal device in which only a small gate voltage would be 
enough to pass through two high/low sequential switching states. 
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